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Starting with a scalar field in a thermal bath using the one loop quantum corrections potential, 
we rewrite the Klein-Gordon equation in its thermodynamical representation and study the phase 
transition of this scalar field due to the U(l) symmetry breaking. We find the corresponding 
thermodynamic and viscosity expressions and propose that these terms can be measured in the 
laboratory in order to know if the phase transition of a Bose-Einstein Condensate can be explained 
in terms of quantum field theory in a simple and elegant way. 
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Introduction.- Symmetry breaking (SB) is one of the 
most essential concepts in particle theory and has been 
extensively used in the study of the behavior of particle 
interactions in many theories [l], 0] ■ 

Phase transitions are changes of state, related with 
changes of symmetries of the system. The analisys of SB 
mechanisms have turn out to be very helpful in the study 
of phenomena associated to phase transitions in almost 
all areas of physics. Bose-Einstein Condensation (BEC) 
is one topic of interest that uses in an extensive way the 
SB mechanisms [3|. The increasing studies in the corre- 
spondence between SB in condensed matter physics and 
cosmology are one of the motivations of this research [4j . 
Although phase transitions happen in condensed mat- 
ter systems and the Universe, their status is not simi- 
lar in both cases. For condensed matter, more transi- 
tions are known, for which our knowledge is better 043 ■ 
Nevertheless, there have been very few experiments con- 
nected to them in other directions likecosmology. The 
fact that condensed matter might permit for experiments 
that open the possibility of examining methods applied 
or used in cosmology and other areas of physics is one 
the most interesting aspects of such research. 

The similarity between transitions in the Universe and 
condensed matter is exciting, but for them to be conclud- 
ing we have the need of studying how transitions take 
place in condensed matter, so that we can see if the same 
assumptions about the Universe take us to understand 
it better. This might be done by looking for character- 
istics of transitions from which we can understand the 
nature of SB. The theoretical expectation that at high 
temperatures phase transitions might be possible in our 
Universe, associated with the breaking of symmetries, is 
of special concern for cosmology. The cosmic field sce- 
narios follows the concept from particle physics and con- 
densed matter physics that as soon as the temperature of 
a system decreases, then a field can develop a nearly def- 
inite classical value that can fix nearly all of its physical 
properties. These field energy concentrations can then 
maybe gravitationally deviate an initially uniform mat- 



ter distribution, possibly seeding structure formation or 
giving other interesting results for the evolution of the 
Universe HQ. 

The results from finite temperature quantum field the- 
ory raise important challenges about their possible man- 
ifestation in condensed matter systems, which can be 
drawn by non-relativistic or relativistic scalar field (SF) 
theories in the framework of the Ginzburg-Landau the- 
ory of condensed systems. Here, our aim is to study a 
relativistic model made up of a real SF together with 
the possibility that this SF dark matter might undergo 
a phase transition as the temperature of the Universe is 
lowered. To analyze this SB we take a model having an 
U(l) symmetry, considering the easiest case of a double- 
well interacting potential for a real SF $(x,t) that goes 
as 
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Here one important idea, to which we shall refer, is that 
of identifying the order parameter ^ which character- 
izes the phase transition with the value of the real scalar 
quantum field $. 

Theoretical framework.- In a expanding universe, the 
SF is governed by the Klein-Gordon (KG) equation given 
by 



□$ + V,$$ <S> + 2F,$ = 0, 



(2) 



where we have added an external interaction potential <fi 
and the D'Alambertian operator which for an expanding 
Universe is defined as 
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being a the scale factor and H = a/ a the Hubble param- 
eter. 

Here we are interested in the properties and behav- 
ior of the SF in a thermal bath, and in the interactions 
with other particles. In that case, the SF potential (p} 
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extended to one loop is Tol 
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cf/h , ks is Boltzmann's 
constant and T is the temperature of the thermal bath, 
this result includes both quantum and thermal contri- 
butions. When the temperature T is high enough, the 
minimum of the potential ([3]) is $ = 0. At this point, 
the SF density is equal to = l/2$ 2 + V{$) ~ T T 4 
and it behaves as a radiation fluid. Now, we suppose that 
the temperature is sufficiently small so that the interac- 
tion between the SF and the rest of matter decouples. 
After this moment, the term TqT is not longer impor- 
tant, thus, for sufficiently low T the radiation term can 
be dropped out. 

The critical temperature where the minimum of the 
potential $ = becomes a maximum is 
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Now for the SF we perform the transformation 
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where k is the scale of the system to be determined by 
the experiment and m is the SF mass in natural units, 
from now on we will take c = 1 . 

In terms of function the KG equation @ reads, 

3 h 2 
-ih{^ + + — - 9A|*| 2 *) 

2 2 m 

+ p*4> - A-^fc|T 2 * = 0, (4) 

where we have used the notation '= d/dt, \ty\ 2 = vp^/* = 
p, and p, = m — \h 2 k 2 B T 2 /Am. Q is KG's equation in 
terms of the function ^ and temperature T. From Q we 
can notice that in a non-expanding Universe, i.e., H = 0, 
when T = and in the non-relativistic limit, □ — > —V 2 , 
eq. Q becomes the Schrodinger equation with an extra 
term \*S>\ 2 . In this limit this is the Gross-Pitaevskii equa- 
tion for Bose-Einstein Condensates (BEC), i.e, an ap- 
proximate equation for the mean- field order parameter. 
The static limit of equation Q is the Ginzburg-Landau 
equation (l2| . 

In what follows we transform the relativistic Gross- 
Pitaevskii like equation Q into its analogous hydrody- 
namical version 
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function "J in its corresponding norm p and a phase S as 
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Here we will interpret p = p/Mx as the rate between 
the number density of particles in the condensed state, 



p = mn — mN /a 3 , being N the number of particles in 
condensed state and Mt the total mass of the particles 
in the system, S = S(x, t) will be a function related to 
the velocity field, as we will see later on, being both of 
them functions of time and position. 

So, from this interpretation we have that when the 
KG equation oscillates around the $ = minimum, the 
number of particles in the ground state is zero, p = = p. 
Below the critical temperature T c close to the second 
minimum, § 2 min = k B (T 2 — T 2 )/4, the density oscillates 
around p = n 2 k 2 B (T 2 - T 2 ) /4. 

Then, from ((4]) and ([5]), separating real and imaginary 
parts we have 
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where w = —9h 2 A/2m 2 . Now, taking the gradient of © 
then dividing by a and using the definition 



v= VS* 
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we obtain, 
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pF^ - - VjJ + pF Q - -pVa, 
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where = — V(/i</>)/ma is the force associated to the 
external potential <j>, p = wp 2 /2 can be seen as the pres- 
sure of the SF gas, Fq = — I/oVJ/q is the quantum force 
associated to the quantum potential [15| . 



Uq 



v 2 Vp 



2m 2 a 2 \ ^fp 
and pVcr/a is defined as 



1 



-pVa = Hvp 1 



1 



CV(lnp) 



-5' 



h 2 p 

Am 2 a 



1 A 



4 m 2 a 



k 2 B pTVT 




(8) 



3 



where the coefficient £ is given by 
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and the term V(lnp) can be written as 
V(lnp) = 



-V(V-v) - ~V[(Vlnp) • v] 
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System (|7|) is completely equivalent to equation 

Observe that from the uncertainty relation mc 2 At ~ %, 
we have that h/(mc 2 )S will only be important when the 
system is relativistic. In the non-relativistic limit, system 
© reads 



'p + -V • (pv) + 3Hp = 0, 
a 

pv + p-(v- V)v = pF - -Wp + pF Q - -pVcr. 
a a a 
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where the V (In p) term now reads 
V(lnp) 



-V(V • v) - -Vf(Vlnp) • vl. 
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and for the non-relativistic limit we have neglected the 
gradient of the term j^t^i ffV Here we interpret the 
function Vcr as a viscosity expression, it contains terms 
which are gradients of the temperature, of the divergence 
of the velocity and of the density, plus a term due to the 
expansion of the universe. 

In what follows we will derive the thermodynamic 
equations from the KG equation. For doing so, we mul- 
tiply (|9b[) by v and divide it by p, with this we obtain, 
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with V = V/a. In the same way, we can derive a conser- 
vation equation for the external potential <j) = pej), given 
by the relationship 

{p$i=pi>+u (io) 

Using the continuity equation (|7al) into (fTU)) we obtain, 

(p4>) + V ■ (pv<j>) + 3Hp4> = -pv • F + p4>. 

Observe that we can generate a continuity equation even 
for Uq in the same way as before, obtaining 

(pU Q J+ V • (pvU Q ) + 3H P U Q = -pv • F Q + pU Q (11) 

Nevertheless, this procedure is not possible for a because 
in general we do not know it explicitly, we only know its 
gradient, only in some cases it is possible to integrate it. 
However, the continuity equation for a would be given 
by 

(per)' + V • (pita) + 3Hpa = pv ■ Vcr + per. 

Now observe that the quantum potential Uq in (fTTj) 
also fulfills the following relationship 

pU Q + V • (pv p ) = 0, 

which follows by direct calculation, and where we have 
defined the velocity density v p by 



v p 4m 2 
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which can be interpreted as the velocity flux of the po- 
tential Uq. Using the continuity equation (|9a|) . this last 
relation can be written as 

(pU Q y + V • [pU Q w + 3 P ) + 3HpU Q = (12) 

where we have defined the quantum density flux 

J p = pv p . 



Equation (|12j) is another expression for the continuity 
equation of the quantum potential Uq. 

On the other hand, the total energy of the system is 
given by 
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where u is the inner energy of the system. By definition, 
the total energy e must be conserved, that means that e 
must fulfill a continuity equations given by 

(pe)'+ V • J e + 3Hpe + p(4> + Uq + cr) +pV • v = 0, (14) 

being J e the total energy current, given by the total en- 
ergy flux and the heat flux, 3 q , 

J e = pew + J q + pv. 
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We now substitute (|13|) into (1141) . in the non-relativistic 
limit we obtain, 

(pu)'+V -(pvu+Jq+pv+J p )+p((j)+a)+3Hpu = —pV-v 

(15) 

In order to find the thermodynamical quantities of the 
system in equilibrium (when p is constant on a volume 
V), we restrict the system to the regime where the grav- 
itational potential is nearly constant in time and where 
the viscosity coefficient fulfills the relation a — 0. With 
this conditions at hand for (fT5l) we have: 



(pu)'+V -{pwu + iq+pw + S^+iHpu = -pV-v (16) 

Multiplying this last equation by a 3 and after integrating 
the resulting expression on a close region, we obtain 

^ J pua 3 dV + <j) a 3 (J q + pv) ■ ndS + <j> a 3 J p ■ ndS 



= -P 



dt 



i 3 dV. 



Equation (JT6J) is the continuity equation for the inter- 
nal energy of the system and as usual, from here we can 
derive a generalized version of the first law of thermody- 
namics for the KG equation or a BEC, this reads 



dU = dQ + AA Q - pdV 



(17) 



where as always, U — J pu a 3 dV is the internal energy of 
the system, Q is the heat, and 



dt 
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<j> a 3 /5(V In p) ■ ndS = - <j> a 3 pv p ■ ndS 



is the corresponding quantum heat flux due to the quan- 
tum nature of the KG equation. The second term on the 
right hand side of equation (TT71) would make the crucial 
difference between a classical and a quantum first law of 
thermodynamics. 

Phase transition.- From hereafter we study the transi- 
tion between the <J> = state to the low energy one with 
T <T C close to the minimum 3> mm = k B \/T 2 - T 2 /2. 

We have already seen that during the epoch when 
T >> T c there are not scalar particles in the ground 
state. We will suppose that the scalar particles decouple 
from the rest of the matter at some moment, such that 
the total density here on remains constant. With the 
space-time expansion the temperature decreases, below 
the critical one T < T c , close to the local minimum the 
density oscillates around the value p = k 2 B n 2 (T 2 — T 2 )/4, 
being ksnT c ~ 1. The function S in (JSJ has a sim- 
ple expression, S — Sot, with sq << rac/h in the non- 
relativistic case. This implies that the velocity Vo = 0, 
also, if there does not exist an external force in the sys- 
tem then, = 0. If the temperature gradients are very 
small, in this limit the density is p — poa~ 3 , thus the 
temperature of the BEC decays as 
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If there is no space-time expansion, the viscosity of the 
BEC might in fact contain the whole information of the 
phase transition. 

Finally, to illustrate the previous exposition, we give 
the following example. Suppose that in the system there 
are only condensed and excited particles of the same 
specie. Thus p = N /(N ex + No), where No is the num- 
ber of condensed particles and N ex the number of exited 
particles. Combining the previous equations we obtain, 
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being N = Nq + N ex . Observe that only a fraction 
N/(kgK 2 T 2 ) of the scalar particles reaches the ground 
state at T — 0, this value can only be determined experi- 
mentally and fits the value of the scale k. The main idea 
we want to point out here is that these phenomena might 
be equivalent for a BEC on earth as for the cosmos, and 
this might follow the previous equations exactly, so this 
function might be tested in the laboratory. If confirmed, 
the phase transition of a BEC can be explained using the 
quantum field theory in a very simple way. We expect 
that all of these results can be checked in an experiment, 
specially the behavior of the density expressions during 
the phase transition of a BEC. 
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